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Abstract. We study families of time-frequency localization operators and de- 
rive a new characterization of modulation spaces. This characterization relates 
the size of the localization operators to the global time-frequency distribution. 
As a by-product, we obtain a new proof for the existence of multi- window Gabor 
frames and extend the structure theory of Gabor frames. 



1. Introduction 

A time-frequency representation transforms a function / on M. d into a function 
on the time-frequency space IR d x M d . The goal is to obtain a description of / 
that is local both in time and in frequency [51 [20] . The standard time-frequency 
representations, such as the short-time Fourier transform and its various modifi- 
cations known as Wigner distribution, radar ambiguity function, Gabor transform, 
all encode time-frequency information. However, the pointwise interpretation of 
such a time-frequency representation meets difficulties because, by the uncertainty 
principle, a small region in the time-frequency plane does not possess a physi- 
cal meaning. Therefore the question arises in which sense the short-time Fourier 
transform describes the local properties of a function and its Fourier transform. 

Following Daubechies [10], we use time-frequency localization operators to give 
meaning to the local time- frequency content. By investigating a whole family of 
localization operators and glueing together the local pieces, we are able to charac- 
terize the global time-frequency distribution of a function. In more technical terms, 
our main result provides a new characterization of modulation spaces. 

We define the short-time Fourier transform (STFT) of a function / e L 2 (R d ) 
with respect to a window function <p £ L 2 (IR d ) as 



(1) V v f(x, uj)= f(t)<f(t - x)e- 2niuJ - t dt, for all z = (x, u) G R 2d . 

The STFT V tp f{z) is a measure of the time-frequency content near the point z 
in the time-frequency plane M. 2d . However, the STFT cannot be supported on a set 
of finite measure by results in [SHIED! SB]- This fact complicates the interpretation 
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of local information obtained from the STFT. In particular, it is impossible to 
construct a projection operator that satisfies V^P^f) = xn • V<pf- As a remedy 
one resorts to the following definition of localization operators. 

We denote translation operators by T x f{t) = fit — x) and time-frequency shifts 
by ir(z)f(t) = e 2niuj4 f(t - x) for x,ui,t G R d . Fix a non-zero function cp G L 2 (R 2d ) 
(a so-called window function) and a symbol a G L x (]R 2a! ). Then the time- frequency 
localization operator H a acting on a function / is defined as 



The integral is defined strongly on many function spaces, in particular on L 2 (R d ). 
A useful alternative definition of H a is the weak definition 



This definition can be easily extended to distributional symbols a G iS'(1R m ). The 
subtleties of the definition and boundedness properties between various spaces have 
been investigated in many papers, see [H [37J EH] for a sample of results. 

If a is non-negative and has compact support in Q C ]R d , then H a f can be 
interpreted as the part of / that lives essentially on Q in the time-frequency plane, 
and so H a may be taken as a substitute for the non-existing projection onto the 
region Q in the time-frequency plane. 

In this paper we investigate the behavior of an entire collection of localization 
operators. Namely, given a lattice A C ~R 2d of the time-frequency plane, we consider 
the collection of operators {H Txa : A G A} and the mapping / — ► {H Txa f}. If the 
supports of T\o cover IR 2d , then {HT x af, A G A} should contain enough information 
to recover / from its local components. In particular, the set {HT xf7 f : A G A} 
should carry the complete information about the global time-frequency properties 
of /. We make this intuition precise and derive a new characterization of modulation 
spaces from it. Similar to Besov spaces, modulation spaces are smoothness spaces, 
but the smoothness is measured by means of time-frequency distribution rather 
than differences and derivatives. Here, we establish a correspondence between 
the behavior of the sequence ||ifr A o-/||2, A G A, and the membership of / in a 
modulation space. 

As a special case of our main theorem we formulate the following result. 

Theorem 1. Fix a non-zero function if in the Schwartz space «S(R d ) and a weight 
function m on M. 2d that satisfies m(zi + Z2) < C(l+ |^i|) N m(z2) for some constants 
C,N > and all z\ : z% G M. 2d . Then a tempered distribution f satisfies 




(2) 




(3) 




if and only if 



(4) 




AeA 
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The expression in ([3]) is just the norm of / in the modulation space M^(R d ). Our 
main result shows that the expression in (using the time-frequency components 
of /) is an equivalent norm on the modulation space M^(M d ). 

In pseudodifferential calculus one often defines spaces by conditions on their time- 
frequency components. For instance, Bony, Chemin, and Lerner [31 H] introduced a 
Sobolev-type space H{m) by using Weyl operators instead of localization operators. 
For the (extremely simplified) case of a constant Euclidean metric on the time- 
frequency plane, a distribution / belongs to H(m), whenever for some test function 

on R 2d 



is finite, where a w is the Weyl operator corresponding to the symbol a. The only 
difference between (jSJ) and (jlj) is the use of Weyl calculus instead of time-frequency 
localization operators and a continuous definition instead of a discrete one. It was 
understood only recently that H{m) coincides with the modulation space M^ l (R d ) 
and that ((Sj) is an equivalent norm on M^(IR d ) [26J. Thus Theorem [1] can be 
interpreted as an extension of [3] to L p -like spaces. 

Let us also mention that in the language of [36], the operators {Ht x( t, A G A} 
form a g-frame for L 2 (IR d ). Our construction seems to be one of the few non-trivial 
examples of ^-frames that are not frames. 

In this paper we prove the norm equivalence of Theorem [1] for a large class of 
modulation spaces and arbitrary time-frequency lattices. For a rather restricted 
class of lattices, namely lattices with integer oversampling, an analogous result 
was derived in [12] for unweighted modulation spaces. The main arguments for 
the integer lattice were based on Zak transform methods and interpolation. For a 
general lattice, these methods are no longer available, and we have to develop a 
completely new approach to some of the key arguments. 

As a by-product of the new techniques we have found several results of indepen- 
dent interest. 

• We formulate several structural results and characterizations of Gabor frames 
for multi-window Gabor frames. 

• We prove a finite intersection property for time-frequency invariant sub- 
spaces of the distribution space M°°(lR d ). This property resembles the 
finite intersection property that characterizes compact sets. 

• We give a new, independent proof for the existence of multi-window Gabor 
frames with well-localized windows. Previous proofs were based on coorbit 
theory [15] and the theory of projective modules [33]. Our proof provides 
additional insight how the windows can be chosen. 

• We derive precise estimates for the localization of the eigenfunctions of a 
localization operator. 

This paper is organized as follows. In Section 2 we recall necessary facts from 
time-frequency analysis. On the one had, we introduce modulation spaces and 
explain their characterization by means of multi-window Gabor frames. On the 
other hand, we state and prove several properties of localization operators. In 
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Section [31 we formulate and prove our main result (Theorem [8]). In Section [37^1 we 
analyze some of the consequences of Theorem [8] and its proof. In the appendix we 
collect and sketch the proofs of some of the structural results on Gabor frames. 

2. Time-Frequency Analysis of Functions and Operators 

2.1. Modulation Spaces. Modulation spaces are a class of function spaces associ- 
ated to the short-time Fourier transform ([1]). Note that for a suitable test function 
ip, the short-time Fourier transform can be extended to distribution spaces by 
duality and V ¥> f(z) = (f,ir(z)<p). 

For the standard definition of modulation spaces, we fix a non-zero "window 
function" g G 5(IR d ) and consider moderate weight functions m of polynomial 
growth, i.e., m satisfies m{z\ + z 2 ) < C(l + |2; 1 |) s m(z 2 ), zi,z 2 G R 2d for some 
C, s > 0. Given a moderate weight m and 1 < p, q < oo, the modulation space 
M^R ') is defined as the space of all tempered distributions / G iS'(R d ) with 
V g f G L™(R 2d ), with norm 

(6) II/IIm&«(r*) = WVgfWi,™^) ■ 

If p = q, we write M^(R d ). 

For weight functions of faster growth we have to resort to different spaces of 
test functions and distributions. Let g(t) = e™ 71 "*'* be the Gaussian window and 
Ti.Q = span {ir(z)g : z G M. 2d } be the linear space of all finite linear combinations 
of time-frequency shifts of the Gaussian. Let v be a submultiplicative even weight 
function on M? d and m be a //-moderate function; this means that v{z\ + z-i) < 
v{z\)v{z-i), u(z) = v(—z) and m(zi + z 2 ) < iy(zi)m(z 2 ) for all z,zi,z 2 G M. 2d . For 
1 < p,q < oo the modulation space M^ 9 (R ) is then defined as the closure of Ho 
in the norm ||/||]vr& 9 (R d ) as in (jBj). If p = oo or q = oo, we take a weak-*-closure 
of TCo. These general modulation spaces possess the following properties. Assume 
that m is //-moderate and 1 < p, q < oo, then 

(7) Ml(R d ) C M p m q {R d ) C M^ u (R d ) = Ml(R d )* . 
Further, if cp G Mj(l rf ), then 

(8) \KfL™~\\vj\\ h ™ = \\f\\ M ™. 

thus different windows in M^(M d ) yield equivalent norms on M^ 3 . 

The embedding (JTj) says that M^(R d ) may serve as a space of test functions 
and M.'^ i/ (R d ) as a space of distributions for all modulation spaces M^; 9 with a 
//-moderate weight m. 

If v s (z) = (1 + \z\) s , s > and m is z/ s -moderate, then we have 

S(R d ) C M^ s (R d ) C Mf^(M a! ) C M^ ! , s (R a! ) C ^'(M ') , 

in agreement with the standard definition, but for v{£) = e a ' 2 ' 6 with a > and 
< 6 < 1 we have 

M^(M d ) C C S'(R d ) C Mw 1/ (M li ) . 
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In the sequel we will start with a submultplicative weight v and take M^ !/ (R a! ) as 
the appropriate distribution space. Our results hold for arbitrary submultiplicative 
weights v. 

For the detailed theory of modulation spaces we refer to [21] Ch. 11-13], for a 
discussion of weights and possible distribution spaces see [23] . 

Sequence space norms. Recall that a time-frequency lattice A is a discrete 
subgroup of IR 2a! of the form A = AH? d for some invertible real-valued 2d x 2d- 
matrix A. 

Given a lattice A C R 2d with relatively compact fundamental domain Q, the dis- 
crete space i^(A) consists of all sequences a = (cia) AeA for which the norm 

(9) ||a|||M = ||X)I aIxa+oIIlSi« 

AeA 

is finite. If A = aL d x 6Z d , then this definition reduces to the usual mixed-norm 
space £^ 9 (Z M ) with norm 

i/p\ 1/9 



^2 ( \ a kn\ p m(ak, bn) 1 



n€Z d fceZ d 

As a technical tool we will need amalgam spaces (in one place only). A mea- 
surable function F on M. 2d belongs to the (Wiener) amalgam space W(L™), if the 
sequence of local suprema 

a kn = essswp XiW £ [0>1] d\F(x + k,u + n)\ = \\F ■ T^xlU 

belongs to £™(Z 2d ). The norm on W(L™) is ||F|| w(l p ;9) = ||a||^. See [27] for an 
introductory article. We need their behavior under convolution and their properties 
under sampling. 

(a) Convolution in Wiener amalgam spaces: Let 1 < p, q < 00 and let m be a 
//-moderate weight. Then 

( 10 ) \\F * < CII^IIwcl^II^IIli- 

(b) Sampling in Wiener amalgam spaces: For F G W(L™) the following sam- 
pling property holds: 

(11) ll-^UIIc 9 - Ca ll-^ll W(L^ 9 )- 

These statements are proved in [27] or [211 Prop. 11.1.4., Thm. 11.1.5.]. 

2.2. Gabor frames. Gabor frames are closely linked to modulation spaces. They 
constitute "basis-like" sets for modulation spaces and are used to characterize the 
membership in a modulation space by the magnitude of coefficients in the corre- 
sponding series expansion. 

For a given lattice A C M. 2d and a window function ip G L 2 (IR d ), let Q(ip,A) 
denote the set of functions {ii(X)ip : A G A} in L 2 (M d ). The operator 



AeA 
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is the frame operator corresponding to Q(<p,A). If is bounded and invertible 
on L 2 (IR d ), then Q(ip,A) is called a Gabor frame for L 2 (R d ). This property is 
equivalent to the existence of two constants A, B > such that 

(12) A\\f\\l<J2\(f^W9)\ 2 =(S v fJ)<B\\f\\l forall feV(R d ). 

AeA 

Using several windows ip = ((fi, . . . <p n ), we say that the union (Jj=i S((fj,A) is 
a multi-window Gabor frame, if the associated frame operator given by 

n n 

(13) S v f = 5252(fM>)<Pi)*W<Pj = XX/ 

3=1 AeA 3=1 

is invertible on L 2 (]R d ). The frame operator can be expressed as the composition 
of the analysis operator C^a defined by 

CV,a(/)(A,j) = (/,7t(A)^->, AeA,i = l,...,n. 

and the synthesis operator -D^a defined by D^ A (c) = J2\eA Y^j=i ^j^W^j- Then 

2.3. Characterization of Modulation Spaces with Gabor Frames. The fol- 
lowing characterization of modulation spaces by means of multi-window Gabor 
frames is a central result in time-frequency analysis and useful in many applica- 
tions. It is crucial for the proof of our main theorem (Theorem |H]). 

Theorem 2. Let v be a submultiplicative weight on M? d satisfying the condition 
lim n ^ 00 v{nz) l l n = 1 for all z G M 2d and let m be a v-moderate weight and 1 < 
p, q < oo. Assume further that UjLi G(<fj, A) is a multi-window Gabor frame and 
that ifj G Ml(R d ) for j = 1, ... , n. 

(i) A distribution f belongs to M.^ n {R d ), if and only ifC Vj f G t v m for j = 1, . . . ,n. 
In this case there exist constants A,B>0, such that, for all f G M^(M d ), 

A\\f\U m < (E (X l(/,vr(A)^,)| 2 ) P/2 m(Ar) 1/P < B\\f\\^. 

AeA j=l 

(ii) Assume in addition that A = aZ d x VL d is a separable lattice. Then a 
distribution f belongs to M^ 9 (IR d ) if and only if each sequence C Vj f(ak,bl) = 
(f,7r(ak,bl)(pj) belongs to £% l q (l? d ). In this case there exist constants A and B 
depending on p, q, m such that, for all f G M^ 9 

(14) 

A\\f\\ MV <E(E (X K/,vr(aA;,6/)^)| 2 ) P/2 m(aA;,6/r) 9/P ) 1/9 < B\\f\\ M ™ . 
zgz fcez 3=1 

^ Hi) Let A C M 2d be an arbitrary lattice and Q be a relatively compact funda- 
mental domain of A. Then a distribution f belongs to M^ 9 (R ) ; if and only if 

( V^ 2 
the function X\eA ( Sj=i I (/> 7r (^)v 9 j) | 2 ) Xa+q belongs to L^ 9 (IR 2d ). In t/ws case 
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there exist constants A,B>0, such that, for all f G M.™ (M. 2d ) , 

n 1 /2 

A||/|| M -< ^(El^^)^)! 2 ) xa+qIIl-^^II/IIm-. 

AeA j=i 

Note that (ii) follows from (iii), since for Q = [0, a] d x [0, b] d the norm equivalence 

II J2k,iez^a k iX(ak,bi)+Q\\^ q x IMIft 9 nolds • 

Theorem [2] has a long history. It extends the basic characterizations of modu- 
lation spaces by Gabor frames to multi-window Gabor frames. For Gabor frames 
with a single window and lattices of the form A = aZ d x b7h d with ah G Q Theo- 
rem [2] was proved in [16j. For general lattices it follows from the main result in [23] 
and the techniques in [16]. See also the discussion in [211 Ch. 13]. The proofs for 
multi-window Gabor frames require only few modifications, we therefore postpone 
a discussion to the appendix. 

2.4. A New Characterization of Mult i- Window Gabor Frames. The proof 
of our main statement relies on a characterization of multi-window Gabor frames 
without using inequalities. The following lemma is a generalization of [22J from 
Gabor frames to multi-window Gabor frames. 

Lemma 3. Assume that ipj G M 1 (M d ) for j = 1, . . . ,n. Then the following prop- 
erties are equivalent. 

(i) \J^ = iG{fj,A.) is a multi-window Gabor frame for L 2 (R d ) . 

(ii) The analysis operator C V) a is one-to-one from M°°(IR d ) to £°°(A, C n ). 

The idea of the proof will be given in Appendix [Aj where we will also list many 
more equivalent conditions. 

2.5. Properties of Localization Operators. We next recall some elementary 
properties of the localization operators Ht x u- Time-frequency localization opera- 
tors have been introduced and studied by Daubechies [11] [TO] and Ramanathan 
and Topiwala [53], and are also called STFT multipliers, time-frequency Toeplitz 
operators, Wick operators, time-frequency filters, etc. They are a popular tool 
in signal analysis for time- frequency filtering or nonstationary filtering [321 E5] , in 
quantization procedures in physics [1], or in the approximation of pseudodifferen- 
tial operators [9], [3Tj . For a detailed account of the early theory we refer to Wong's 
book [39], for a study of boundedness and Schatten class properties to [71 151 [T8| [37] . 

Lemma 4 (Intertwining property). // a G L°°(M M ) ; ip G L 2 (R d ), and A G A, then 

ix{\)H a <K{\y = H T ^. 

The proof consists of a simple calculation, see [T2l Lemma 2.6]. 

For estimates of the STFT of H a f we introduce the formal adjoint of V v , namely 

V;F= [ F{z)7t{z)cpdz, 

JR 2d 
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which maps functions on R 2d to functions or distributions on R d . With this notation 
we can write the localization operator H a as 

H a f = V;(aVj). 

The STFT of V*F satisfies a fundamental pointwise estimate [21, Proposition 
11.3.2.]: 

(15) \Vv(V;F)(z)\ < (\VM * |F|)(z) Vz G R 2d . 
We note that for F = cV v f this estimate becomes 

(16) K(H a f)(z)\ = K(v;(avj))(z)\ < (\vM*(*Kf\))(z). 

Thus the short-time Fourier transform of H a is a so-called product-convolution 
operator. The standard boundedness results for localization operators can be easily 
deduced from the well established results for product convolution operators [BJ. 

Estimate (fTB"|) is quite useful for the derivation of norm estimates. In the 
following, we fix a non-negative symbol a and investigate the set of operators 
{Hr x a : A G A}. To simplify notation we will write H\ instead of Ht xU , and 
sometimes Hq = H a by some abuse of notation. 

Lemma 5. (%) Assume that a G L^R 2 ^), o > and that ip G L 2 (R d ). Then each 
H\, A G A, is a positive trace-class operator. 

(ii) If, in addition, <p G M^(R d ) and a G L^R 2 ^), then each H\ is bounded 
from M°°(R a! ) into M^(R d ). In particular, all eigenfunctions (pj of H a belong to 
Ml(R d ). 

(Hi) Furthermore, if (p G M*(R d ) and o G L^(R 2d ), then each H\ is bounded 
from ~M.y v (R d ) into L 2 (R d ). 

Proof. Statement (i) is well-known, see, e.g., [2| [T71 [39] . 
To show (ii), we use (TJBD to obtain, for / G M°°(R' i ), 

WUl = llV^/fc = ||V v V;(aV^)|| Li 

(17) < |||v^|*kv v /||| L1 

< ||v^|| L j lkv v /|| L j , 

where we have used Young's inequality. Since <p G M^(R d ) if and only if V^ip G 
L^(R M ) by [2B Prop. 12.1.2], we find that 

||#W||mi < ||V v ¥>||li lkllLil|V^/||L» < C||«7|| L i||/||m« , 

and thus H a is bounded from M°°(R d ) to M^(R d ). 

The proof of (hi) is similar. Again, we apply (1161) to obtain for / G M^^(R a! ): 

\\H*f\\v < \\\VM*\<rV<rf\\\u ^ W V M\v \W V<pf\\v- , 
Hence, the result follows from 

J R 2d V\Z) 

□ 
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The spectral theorem for compact self adjoint operators provides the following 
spectral representation of H x . 

Corollary 6. Assume <p G M^(IR d ) and a G L^IR 2 ^). Then there exists a positive 
sequence of eigenvalues c = (cj) G i 1 and an orthonormal system of eigenf unctions 
(fj G M^(IR d ), such that 

oo 

(18) H a f = J2 c ^^r 

3=1 

It follows that 

oo 

(19) H x f = H Txa f = 7r(A)^7r(A)V = £ Cj (f, 7r(A)^)vr(A)^, 

3=1 

and {n(\)tpj,j G N} is an orthonormal system of eigenfunctions of H\. 

A priori, the spectral representation of H\ holds only for / G L 2 (R d ). The next 
corollary extends the spectral representation to all of M.y u (JBL d ). 

Corollary 7. The expansion for Hxf given in ( |19p is well-defined on M^ i/ (M d ) 
and converges to H x f in L 2 for all f G M.y u (E. d ). 

Proof. Without loss of generality, we assume A = and set H = H a . Since 
Hf G L 2 (R d ) for every / G M^(R d ) by Lemma E3(iii) , we can expand Hf with 
respect to the orthonormal system of eigenfunctions of H and obtain that 

oo 

(20) Hf = J2(Hf,^+r 

3=1 

for some r G L 2 (R d ) in the orthogonal complement of span {<£>.,• : j G N}. As 
H is self-adjoint on L 2 (R d ), we also have (Hf,ipj) = (f,Hipj) = Cj(f,ipj), and 
consequently 

oo 

(21) Hf = ^2c s (f i ip j )i Ps + r. 

3=1 

We need to show that r = 0. Since r G L 2 (M a! ) is orthogonal to all eigenfunctions 
(fj, we find that (Hf,r) = ||r|||. 

To show r = 0, we first observe that (Hh,r) = for all h G L 2 (M d ) by ([IB"]) . 
Since L 2 (M d ) is w*-dense in M^ i/ (IR d ), we may choose an approximating sequence 

f n G L 2 (M d ) such that f n ^ f G M^ jy (M d ). For instance, /„ may be chosen as 

/»= / XB n {z)V g f{z)'K{z)gdz, 

where B n is the ball with radius n and centered at 0. Furthermore, since f n — > f, 
we obtain in particular that V v f n converges to V v f uniformly on compact sets 
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Theorem 4.1]. Consequently 

= (Hf n ,r) = [ a{z)V v f n {z)V~F{z)dz - / a(z)VJ(z)V^)dz = (Hf,r) = 

This shows that r = and so the series (fl9l) represents Hf for all / G M^ jy (M d ). □ 

3. From Local Information to Global Information 

We first state and prove the main result for the modulation spaces M^(IR d ). The 
generalizations to M^ 9 (M d ) will be discussed later. As always, v denotes a sub mul- 
tiplicative, even weight function on M? d satisfying the condition lim^oo v{nz) l ^ n = 

1 for all z G R 2d . 

Theorem 8. Let a G Lj(R M ) be a non-negative symbol satisfying the condition 

(22) A < ^T x a < B, a.e. 

AeA 

for two constants A, B > 0. Assume that <p G Mi(M d ). Then for every v-moderate 
•weight m and 1 < p < oo the distribution f G M^(lR d ) belongs to M^(lR d ) ; if and 
only if 

(23) (j]||i7 A /||XAr) 1/P <oo, 

AeA 

and the expression in ( 123]) is an equivalent norm on M^(IR d ). 
Similarly, for p = oo we obtain the norm equivalence 

(24) ||/|| Ms xsup||F A /|| 2 m(A). 

AeA 

The norm equivalence supports the interpretation that H\f carries the local 
time-frequency information about / near A G M. 2d . By combining the local pieces 
H\f, one obtains the global time-frequency information as it is measured by mod- 
ulation space norms. 

The proof of Theorem [8] requires some preparations. We first show that finitely 
many eigenf unctions of Hq = V^crV^ generate a multi-window Gabor frame for 
L 2 (M (i ). With this crucial step in place, Theorem [8] can then be deduced from the 
characterization of modulation spaces by means of Gabor frames. 

3.1. Multi- Window Gabor Frames. 

Lemma 9. Assume that a G L 1 ^) and ^2 X£A T x cr x I, and that <p G M^lf^). 
Let {ifij : j 6 N) be the orthonormal system of eigenfunctions of H . Then there 
exists n G N, such that the finite union Uj=i 0((Pj, A) is a multi-window Gabor 
frame for L 2 (IR rf ). 

An analogous statement was proved and used in [12] for the lattice A = 7? d and 
rational lattices by means of Zak transform methods. In the case of general lattices 
we cannot apply Zak-transform methods. As a substitute, we will use a finite 
intersection property for A-invariant subspaces of M°°. The following statement 
may be of interest in its own right. 
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Lemma 10. Assume that W n is a sequence of w*- closed subspaces in M°°(IR d ) 
such that 

(i) W„ 2 W„+i ^ {0} /or all n EN and 

(ii) W n invariant under all operators 7r(A) /or A G A. 

T^en fl„>i>V^{0}. 

Proof. Let Q be the closure of a relatively compact fundamental domain of A, for 
instance, if A = AZ 2d , then Q = A[0, l] 2d . We first choose a sequence h n G W„ 
with ||/i n ||M°° = su P2giR2d IV^/i^z)! = 1. Then there exists a sequence of points A„ 
in A, such that 

sup |V^(7r(A n )/i„)(z)| = 1 . 

zeQ 

Since W n is invariant under all 7r(A),A G A, the distribution /„ = Ti(\ n )h n is in 

Next we show that the set of restrictions {V^fnlo) ls equicontinuous. We have 
(25) Kf n {z)-VM0\ = \(fn,(ir{z)-7r{0)<p)\ < H/Jm- ■ ||(vr(^) -7r(£)M| M i • 

Since ||/ n ||M°° = H^ArO^-nllM 00 = 1, the equicontinuity follows from the strong 
continuity of time- frequency shifts on M 1 (R a! ). 

We next choose z n G Q with iV^/n^n)! > |- Since the unit ball in M 00 (M <i ) is 
w*-compact, there exists a subsequence f nk that converges to some / G M°°(M d ) in 
the w;*-sense. Furthermore, by compactness of Q, there also exists a subsequence 
Zi of z nk , such that ^ — > z G Q. Hence, by equicontinuity, 

V v /,(z,) - V„/(z) . 

Since \V v fe(ze)\ > 1/2, we conclude that also \V tp f{z) \ > 1/2, and consequently 
/ / 

By construction, fi G W m for every £ > m, hence we obtain / = w*— lim^oo fi G 
W m for all m, because W m is w*-closed. To summarize, we have constructed a non- 
zero / G M°°(lR d ) that is in W m for all m. 

□ 

Proof of Lemma To prove that finitely many eigenfunctions generate a multi- 
window Gabor frame with respect to the lattice A, we assume on the contrary that 
Uj=i A) is not a frame for every n G N. Using Lemma [TU1 and Lemma [21 we 
will derive a contradiction to the assumption that A < X^AeA-^ " < B. 

We use the criterion of Lemma [31 Let (fi n = ((pi, . . . , (p n ) be the vector- valued 
function consisting of the first n eigenfunctions of H , and 

W„ = ker(C VnjA ) = {/ G M°°(M d ) : (/, tt(A)^) = 0, VA G A, j = 1, . . . , n} 

be the kernel of the coefficient operator C Vni A in M°°(]R rf ). 

If Uj=i ^ (Vj , A) is not a frame, then W n is a non-trivial subspace of M°°(IR d ) by 
LemmaO By construction, the W n 's form a nested sequence of w*-closed subspaces 
of M°°(lR d ), and they are also invariant under vr(A), A G A. Thus the assumptions 



12 



MONIKA DORFLER AND KARLHEINZ GROCHENIG 



of Lemma [TU1 are satisfied, and we conclude that f]^Li W n 7^ {0}- This means that 
there exists a non-zero / G M°°(R d ), such that 

(26) (/, n(\)<Pj) = for all A G A and all j G N. 

We now consider H\f . Since H\f G M 1 (IR d ) by Lemma the bracket (H\f, f) is 
well-defined and given by 

(27) (H x fJ}= [ a(z-X)\VJ(z)\ 2 dz. 

On the other hand, the extended spectral representation of Lemma [7] and (j2~6l) 
imply that 

(28) H\f = ^ Ci (/,7r(A)^-}7r(A)^ = 0. 

Consequently (H\f, f) = for all A G A, and | V^,/" (2;) | 2 vanishes on [J XeA supplier. 

According to the crucial assumption (122]) we have YlxeA^xv — ^ > ^ almost 
everywhere, and thus U AeA supp(T\er) = IR 2d . Therefore, (1271) and (1281) imply that 
Vipf — 0, from which / = follows. This is a contradiction to / being a non-zero 
element in fl^Li W n . 

This contradiction shows that there exists an n G N, such that (Jj=i ^(fj^) * s 
a multi-window Gabor frame, and we are done. □ 

Remark 1. Note that for finite-rank operators i?o ; h can be seen directly that the 
finite set of eigenvectors generates a multi-window Gabor frame for A. 

3.2. Proof of Theorem [8]. We are now ready to prove the main theorem. We 
observe that for / G M^(lR d ), H\f G L 2 (R d ) by Lemma Mni) . Thus the terms in 
( T2"3"j) are well-defined. 

First assume that p < 00 and / G M^(IR d ) C M™ /u (R d ). Using the embedding 

M 1 ^) <-> L 2 (R d ) and the estimate (HZ) with z/ = 1, we majorize ||#a/||2 as 
follows: 

(29) \\H x f\\ 2 < CjH x f\\ M1 

= C^C [ \a(z-X)\-\V 9 f(z)\dz 
= ^C(|V^|*a v )(A), 

where u v (z) = a y (—z). Thus ||iT\/||2 is majorized by a sample of \V v f \ * a y . 
To proceed further, we use the fact that V v f G VF(L^) and || V^f\\ w(l4) — 
CoMMI/Hm* for up G Ml(R d ) and / G M^(M d ) by [21, Thm. 12.2.1]. Now 
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the convolution relation ffTOl) and the sampling inequality (fill) imply that 
(30) E WxfWSmW < C V C\\(^ * \V v f\)\ A \\ p £Pm 



<C v CC A \\^*\V v f\\\w^ m ) 

< C 9 CC A \\a\\l4VM P w{ ^) ^ ^^AlMiyi/IlM- 
The same argument yields sup AeA ||ifx/|| 2 ?7i(A) < C||/||m^- 

Hence, for 1 < p < oo, the mapping / — > (||iJ A /|| 2 )AeA is bounded from M.^ n (R d ) 
toC(A). 

Conversely, assume that p < oo and 

^||i/ A /||fm(A) p <oo. 

A 

We need to show that / G M^(M d ). Since ||#a/||2 = su P|| s || 2 =i \(H\f,g)\, we have 
the inequality 

\(H x f,g x )\Pm(\y < \\H x f\\ p 2 m(\y < oo 

A A 

for arbitrary sequences g x G L 2 (IR d ) with \\g\W2 = 1- Applying the eigenfunction 
expansion of Corollary El we obtain 

00 p 
( 31 ) E J2 c i(f>< X ^^ X ^>9x) m(\r <J2\\Hxf\\ P MXY < 00. 

A j=l A 

Now fix jo ^ N and set g\ = 7r(A)(^j for A G A. Since the eigenfunctions of H x are 
orthonormal, the sum over j collapses to a single term, and fl3Tl) becomes 

£ |<# A /,2 A )|MA) P = E \c jo {fM^io)\ P ™(X) p < E II^A/||fm(A)' < oo. 

A A A 

The last inequality holds for every j G N. After summing over finitely many j 
and switching to the £ 2 -norm on C n , we obtain the inequality 

1/2 



E ( E 1 </> *w<pj) 1 2 ) m ( A ) p ^ E E 1 (/. i Pm ( A ) ? 

A i=l j=l A 

(32) < *) WHxfWlmiXY < 00. 

i=l ^ A 

We now apply Lemma [9] and choose an n G N, such that Uj=i ^(^-, A) is a multi- 
window Gabor frame for L (R ). Since all (fj are in Mj(R d ), the fundamental 
characterization of modulation spaces (Section 12.31) is valid. Thus Theorem |2]^i) 
implies that / G M p m (R d ). 

If p = 00 and sup AgA || || 2 ""T-(A) < 00, then, by choosing g x as before, we find 



c J0 sup \ (f,ir(X)(p jo )\m(X) < sup \\H x f\\ 2 m(\) < 00 

A A 



for every j . 



M° 
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Arguing as above, Theorem [2] says that 

<C max sup | (/, ir(\)<pj) |m(A) < ( max — ^ sup \\H\f\\2m(\) < oo. 

j=l,...,n x \ j=l,...,n Cj / A 

and / e M^(R M ). 

Combining ( 1301) and (1321) . we have shown that ||/||m^ an d (Eaca II^A/||2 m (^) 
for 1 < p < oo (or sup AeA ||i?A,/||2 m (A) for p = oo) are equivalent norms on 

3.3. Variations of Theorem [5]. In order to formulate our main result for mixed- 
norm spaces and arbitrary lattices, we have to resort to the theory of coorbit spaces, 
as introduced in [131 [H]. in particular, for arbitrary lattices, a sequence (ca)aca is 
in the sequence spaces associated with L^ 9 (R 2d ), if J2\eA c aXa+q is in L^R 2 ^) for 
some fundamental domain Q of A. With this definition, we may give the following 
characterization. 

Theorem 11. Let A be an arbitrary lattice in R 2d and Q be a relatively compact 
fundamental domain Q. Assume the same conditions on o and (p as in Theorem^ 
Then a distribution f G My v (R d ) belongs to M% l q (R d ), 1 < p,q<oo, if and only 

if 

(33) £||#a/|| 2 Xa +0 GL-(R m ), 

AeA 

and || EacA II^a/||2XA+q|| l p, 9 x II/IIm™- 

Proof. The proof is almost identical to the proof of Theorem [SJ The only modifi- 
cations occur in fl30|) . which has to be replaced by 

II \\ H *fhx\+Q\k™ < \\^2\v<pf * °W\x\+q\\-l™ < c\\v v f * ■ 

AeA AeA 
Likewise, in (1321) we replace the weighted L^-norm by the general L^-norm. □ 

For a separable lattice A = aZ d x b7L d the norm in ( 1331) is just the ^f-norm on 
Z with m(k,n) = m(ak,bn). In this case, A = (ka,nb), k,n £ Z d and we may 
write H x f = H k>n f. 



i/p 



Corollary 12. Let A = aL d x b7L d be a separable lattice and assume the same 
conditions on a and tp as in Theorem^ Then a distribution f 6 M^„(l ) belongs 
to M% l q (R d ) for 1 < p, q < oo, if and only if 

( 34 ) ( E(E \\H k ,nf\\ P 2 m(ka,nbr) 9/p ) 1/q < oo , 

n&L d k&L d 

and (134")) defines an equivalent norm on M^' 3 (R d ). The result holds for p = oo or 
q = oo with the usual modifications. 
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3.4. Existence of multi-window Gabor frames and properties of the eigen- 
functions (pj. We finally point out some immediate consequences of our results 
and methods. 

The intermediate results leading to Theorem [8] also imply the existence of multi- 
window Gabor frames for general lattices. 

Theorem 13. Let A be an arbitrary lattice and v a submultiplicative weight on R 2d . 
Then there exist finitely many functions tpj G M^(R d ) ; such that Uj=i G(<Pji A) ^ s 
a multi-window Gabor frame for L 2 (R d ). 

Proof. Choose a G (R 2d ) such that J2 XeA T x a x 1 and fix a window up G (R '). 
For instance, one may choose the characteristic function xq of a (relatively com- 
pact) fundamental domain of A and the Gaussian window (p(t) = e~ nt4 . 

Now consider the localization operator H = V*aV v . According to Lemma E^ii), 
all eigenf unctions tpj of Hq belong to M^(R d ). Lemma [9] states that for some finite 
n G N the set Uj=i ^{Vji ls a multi-window Gabor frame for L 2 (R d ). □ 

The existence of multi-window Gabor frames for general lattices was known be- 
fore. On the one hand, it is an immediate consequence of coorbit theory applied to 
the Heisenberg group. To be more precise, according to [151 Thm. 7] for every lat- 
tice A and every non-zero g G M^(R d ) there exists n G N, such that the set Q(g, -A) 
is a Gabor frame for L 2 (R d ). Using a coset decomposition ^A = U(// + A) for suit- 
able /i G A, one sees that Q(g, -A) = [J Q{ji{n)g, A) is a multi-window Gabor frame 
with all windows Tc(fi)g derived from a single window g. Recently Luef [23] proved 
the existence of multi-window Gabor frames by exploiting a connection between 
Gabor analysis and non-commutative geometry. Our methods provide a third, in- 
dependent proof for this interesting result. 

The construction of multi-window Gabor frames in Proposition [13] yields more 
detailed information about the frame generators, since they are eigenfunctions of 
a localization operator. Intuitively the eigenfunctions corresponding to the largest 
eigenvalues of a localization operator concentrate their energy on the essential 
support of the symbol a of H . For the special case of compactly supported a, this 
intuition is made precise by the following result. 

Proposition 14. Let the non-negative function o G L 1 (R 2d ) be supported in a 
compact set Q in R M with < cr(z) < C a < oo for z G VL. Consider the localiza- 
tion operator given by H a f = V*aV v f with <p G M 1 (R d ) ; ||^|| 2 = 1 and spectral 
representation as in Corollary Then the eigenfunctions (pj of H a satisfy the 
following time-frequency concentration 



(35) 




Equality holds, if and only if a(z)/C a = Xn(z) is the characteristic function of 
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Proof. Using the weak interpretation of H a from (jSJ), we obtain 

\V v <Pj(z)\ 2 dz>-±- [ a(z)\V^(z)\ 2 dz 



□ 



Appendix A. Characterizations of Modulation Spaces and 

Multi- Window Gabor Frames 

In the appendix, we will sketch the proof of Theorem [2] and formulate a series of 
new characterizations of multi-window Gabor frames. These statements generalize 
well-known facts from Gabor analysis and the results about Gabor frames without 
inequalities in [22] . 

For the investigation of multi-window Gabor frames we need the dual concept 
of vector-valued Gabor systems. In this case we consider the Hilbert space TC = 
L 2 (R d ,C") consisting of all vector-valued functions f(t) = (fi(t),...,f n (t)) with 
the inner product 



(36) (f , ^> L 2(Rd,c«) = / fj( t ) l fA t ) dt = X^'^ L2 ( Rd )- 

3=1 J 3=1 

Time-frequency-shifts act coordinate-wise on f. The vector-valued Gabor system 
G((p,A) = {7r(A)<£ : A G A} is a Riesz sequence in L 2 (lR d ,C n ), if there exist 
constants < A, B < oo such that for all finitely supported sequences c, 

(37) A||c|| 2 < || ^2c^7r{X)ip\\l 2(RdjCn) < B\\c\\ 2 2 . 

AeA 

We now proceed to the proof of Theorem [21 The crucial step is to show the 
invertibility of the frame operator on M^(M d ). This step requires a special repre- 
sentation of the frame operator due to Janssen [29] and at its core uses 'Wiener's 
lemma for twisted convolution" [25] . 

For ifj, (f>j in M 1 (lR' i ), j = 1, . . . , n, we denote frame-type operators by 

n n 

s v ,i,f = ^2^2(f, 7r(A)^xA)^ = ■ 

AeA j=l 0=1 

The frame operator of the Gabor system \J^ =1 Q(fj, A.) is S = S^^. We usually 
omit the reference to the lattice A and the windows ipj. 

The volume s(A) of a lattice A = Al? d is defined as the measure of a fundamental 
domain of A and is | det(A)|. The adjoint lattice of A is A° = {// G M. 2d : 7r(A)7r(/x) = 
7r(/i)7r(A) for all A G A}. 

Lemma 15 (Janssen's representation). Assume that tpj,ipj G M 1 (M d ) for all j = 
1, . . . , n. Then the frame type operator associated to Uj=i GiPji-ty an d U?=i ^C0i> A) 
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can be written as 

n 

(38) S„+f = s(A)- 1 I>i. tWW)/ 
mii/i unconditional convergence in the operator norm on L 2 . 

Proof. By Janssen's result [29] the representation holds for a single S^.^ and ( 1381) 
follows by taking a sum. □ 

The canonical dual frame is defined to be 7^ = ^(A)^ -1 ^' Since the frame 
operator S = commutes with time-frequency shifts on A, we obtain the re- 
construction formulas 

n 

/=5- 1 5/ = ^^(/, 7 r(A)^)7r(A) 7j 
AeA j=i 

n 

=S3~ 1 f = ^2^2(f, tt(A)7 j >(A)^ 
AeA 3=1 

As a general principle the localization of a frame is inherited by the dual frame |19j . 
The following statement is a generalization of [25j Thm. 9] to multi-window Gabor 
frames on general lattices. 

Lemma 16. Assume that v is a submultiplicative, even weight on M? d satisfying 
lim n ^ oc v{nz) x l n = 1 for all z G M? d . Assume further that Uj=i ^(v 9 ?; A) is a frame 
for li 2 (M. d ) and that <fj G M^IR 6 '). Then the frame operator S is invertible on 
Ml(R d ) and jj = S^Vj e M^(M d ) for j = 1, . . . , n. 

Proof. Janssen's representation (138]) implies that 

(39) 5 = ^ V5 = S (A)- 1 ^c M 7r( M ), 

AiGA° 

with a coefficient sequence = YTj=i{ { Pj^' K i.t 1 ) { Pj)- The hypothesis ^ G M^(]R d ) 
guarantees that X^eA° l(v 9 j! 7r (A t ) ( ^i)l z/ (A i ) < 00 f° r each j, see [21] Cor. 12. 1.12], and 
therefore the coefficient sequence (c M ) is in ^(A°). Since Uj=i ^(Vj) A) i s a frame, 
the frame operator S VtV is invertible on L 2 (IR rf ). It follows from [25l Theorem 3.1] 
that the inverse frame operator S~ x is again of the form S 1-1 = X^ga° ^ 7r (/ i ) with 
a coefficient sequence d in £l(A°). This representation implies that S 1-1 is bounded 
on M^(lR d ) and that 

(40) INImi = II^VjIImi < cikillMi • 

Therefore the dual windows jj, j — 1, . . . , n are in M*(lR d ) as claimed. □ 

Once the invertibility of the multi-window frame operator on M^(M d ) is estab- 
lished, the proof of Theorem [2] is straight-forward by using the following bounded- 
ness properties of the coefficient operator and -D^a from [2TJ Theorem 12.2.3. 
and 12.3.4.]. If ^ G Ml(R d ) and jj G ' M^(M d ),' then both C VtA and C 7iA 
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are bounded from M^(R d ) into £f^(A,C n ) for 1 < p, q < oo and for every v- 
moderate weight m. Likewise and -D 7 a are bounded from £%?(A, C n ) into 

M^ 3 (M. d ) . For the £^ q {A, C n )-norm we use the Euclidean norm on C n , so that 

ll c ll£fk 9 (A,C") = II EagA (E"=l IcAj'I 2 ) 1 2 Xa+q||l^- 

As a consequence, the reconstruction formula / = D^C^f = -D^a^-va/ holds 
for / G M^(R d ) with the correct norm estimates. The norm equivalence stated in 
Theorem [2] then follows from 

ll/llM^(R d ) = ||-D-y,AC ¥ >,A/||Ma 9 (K d ) - II^,a||opI|£<m/II^ 9 (A,C») 
< 1 1 -D-r, A 1 1 op 1 1 Ccp, A 1 1 op 1 1 / 1 1 M^ 9 (R d ) • 



Next we come to the characterization of multi- window Gabor frames (Lemma [3]) 
and extend the list of equivalent conditions. For the formulation of the dual condi- 
tions on the adjoint lattice A° we need the vector-valued versions of the analysis and 
synthesis operators. For f = . . . , f n ) G M°°(R d , C n ) and <p = (</? x , . . . , (p n ) G 
M 1 (IR d , C n ) the coefficient operator is defined to be C v ^(i)[jj) = ((f, 7r(fx)ip)), /i G 
A°, and the synthesis operators is D^o^c) = X^ga <V 7r u )V- The Gramian op- 
erator G^a = C^aoD^ao is defined on sequences indexed by A°. 

Lemma 17. Assume that ipj G M 1 (IR d ) for j — 1, . . . ,n. The following are equiv- 
alent for the multi-window Gabor system Uj=i S((fj, A): 

(i) [f j=l Q{ip 5 , A) is a frame for L 2 (R d ). 

(ii) Wexler-Raz biorthogonality: There exist jj G M 1 (IR d ),j = l,...,n, such 
that 

n 

(41) s(A)- 1 ^2((pj, 7r(/i)7i) = <W for /i G A°. 

(iii) Ron-Shen duality: Q(tp,A°) is a Riesz sequence in L 2 (R d ,C ?1 ). 

(iv) S VtV> is invertible on M 1 (IR d ). 

( v ) is invertible on M°°(R d ). 

(vi) Stp jip is one-to-one on M.°°(R d ). 

(vii) The analysis operator C<^a : M°°(R a! ) i— > £°°(A,C n ) is one-to-one from 
M°°(R d ) to £°°(A,C n ). 

(viii) The synthesis operator D^ ^ defined on £ X (A, C n ) has dense range in M 1 (IR d ). 

(ix) -D^a is surjective from £ 1 (A,C n ) onto NL 1 (R d ). 

(x) The synthesis operator -D^a defined on £°°(A°) is one-to-one from £°°(A°) 
to M°°(R d ,C n ). 

(xi) The analysis operator C Vt A° defined on M 1 (R d ,C n ) has dense range in 
£\A°) . 

(xii) C^a is surjective from M 1 (R d , C n ) onto£ 1 (A°). 

(xiii) Gtp^o is invertible on £ 1 (A°). 

(xiv) Gtp^o is invertible on £°°(A ). 

(xv) G^a is one-to-one on£ 1 {A°). 
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The equivalence (i) (vi) is claimed in Lemma [3] and is all we need for the 
main results of our paper. 

Proof. The implication (z) =>- (iv) was sketched in Lemma [16j 
(z) <^ (n): Time-frequency shifts on a lattice are linearly independent in the fol- 
lowing sense: if c = (c m ) mG ao g £°° and ^ MgA ° c M 7r(//) = (as an operator from 
M 1 (R d ) to M°°(M d )), then c M = for all // 6 A°, see [22]. Now, if / = 5 ViT / for 
all / G M 1 (IR d ), then by Janssen's representation ( |38l) we have 

n 

/ = ^ ^(^,7r(/i) 7j )7r(/i)/. 

/igA° j=l 

The linear independence of time- frequency shifts implies (141 p . The converse is ob- 
vious. 

(ii) -vv- (iii): Assume first that Ui=i G((fj, ^) ls a multi-window Gabor frame for 
L 2 (R d ). The upper bound in (13 7p follows from the boundedness of the synthesis 
operator on L 2 (IR d ). To show the existence of a lower bound, we apply the 
Wexler-Raz relations. Since Uj=i G(<Pji ls a frame with dual {J™ =1 G('Jj, A.) anc ^ 
7 i G M 1 (M d ) for all j, we have (if,7r(n)j) = Y.]=i( L Pj^ 7r (l J 'hj) = S ( A )<W, and 
A°) and therefore £(7, A°) are biorthogonal systems in L 2 (M d ,C n ). If / = 
X^eA° vOOVj then = s(A)" 1 (f , 7r(/i)7) L a (R d >C n) and 

c = S (A)- 1 C ; ^A°f, 

from which the lower bound in (l3"Tj) follows. 

Conversely, assume that Q(cp,A°) is a Riesz sequence in L 2 (IR d ,C n ). Then 
there exists a biorthogonal basis of the form {7r(/z)7 : \i G A°} contained in 
/C = span((?(<£>, A°)). It can be shown that 7 G M 1 (R d ,C™). The frame prop- 
erty of G{<Pj, A) follows from the Wexler-Raz relations (14TI) . 

With three classical statements (1381) and (ii), (iii) for multi- window Gabor frames 
the remaining equivalences follow exactly as in [22]. □ 
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